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Abstract 

The fundamental biological processes of development of tissues and organs in multicellular or¬ 
ganisms is governed by various signaling molecules, which are called morphogens. It is known 
that spatial and temporal variations in concentration profiles of signaling molecules, which are fre¬ 
quently referred as morphogen gradients, lead to cell differentiation via activating specihc genes in 
a concentration-dependent manner. It is widely accepted that the establishment of the morphogen 
gradients involves multiple biochemical reactions and diffusion processes. One of the critical ele¬ 
ments in the formation of morphogen gradients is a degradation of signaling molecules. We develop 
a new theoretical approach that provides a comprehensive description of the degradation mecha¬ 
nisms. It is based on the idea that the degradation works as an effective potential that drives the 
signaling molecules away from the source region. Utilizing the method of hrst-passage processes, 
the dynamics of the formation of morphogen gradients for various degradation mechanisms is ex¬ 
plicitly evaluated. It is found that the linear degradation leads to a dynamic behavior specified by 
times to form the morphogen gradients that depend linearly on the distance from the source. This 
is because the effective potential due to degradation is quite strong. At the same time, the nonlin¬ 
ear degradation mechanisms yield a quadratic scaling in the morphogen gradients formation times 
since the effective potentials are much weaker. Physical-chemical explanations of these phenomena 
are presented. 


I. INTRODUCTION 


The development of multi-cellular organisms is one of the most important fundamental 
processes in nature [l-3|. The most critical question here is how a small set of genetically 
identical cells in embryos can produce morphologically and functionally different tissues and 
organs in fully developed organisms. The central concept of biological development is that 
the observed complex spatial patterning is a result of action of signaling molecules that are 
also called morphogens [l|-8|. Signaling molecules can produce non-uniform concentration 
prohles, the so-called morphogen gradients, that via complex biochemical networks stimulate 
or suppress specihc genes in embryo cells, depending on the local concentration. In recent 
years, there were multiple experimental and theoretical investigations on how the morphogen 
gradients are created and how they function. This led to several exciting discoveries in the 


held j5-l24|. However, many aspects o 


gradients remain not fully explained 251. 


the mechanisms for formation of the morphogen 


A large variety of approaches to describe the development of morphogen gradient have 
been proposed and discussed 7, 25|. Many of them follow the original idea of Turing that 
the morphogen gradients are resulting from complex reaction-diffusion process 2q. The 
most popular and widely utilized method to explain the formation of sig naling molecules 

1, 3 . In this pic- 


prohles is known as a synthesis-diffusion-degradation (SDD) model 
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ture, the process starts with morphogens being produced at specihc localized regions in the 
embryo, from which they dihuse along the cells. Signaling molecules also can be removed 
from the system after binding to specihc receptors on cells. At large times, this leads to expo¬ 


nential decaying concentration proh 


observed morphogen gradients 
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1 qualitatively agree with many experimentally 


28|. 


It is widely accepted that the process of degradation or removal of signaling molecules 
from the system is critically important for the development of morphogen gradients . This 
allows the formation of the stationary prohles of signaling molecules, ensuring the robustness 
of the genetic information transfer in biological development. But specihc details of how 
the degradation inhuences the formation of morphogen gradients are still not well clarihed. 
There are many counter-intuitive observations that cannot be explained by current theoret¬ 
ical views. In the classical SDD model it is assumed that the degradation is linear, i.e., the 
particle hux leaving the system is proportional to the local concentration of morphogens. 






































It was shown theoretically then that for this model the time to establish a stationary mor- 
phogen gradient at given location, which is also known as a local accumulation time (LAT), 


is a linear function of the distance from the source 


20| . This observation is surprising since 


for the system with unbiased diffusion of particles much more slower quadratic scaling was 
expected 


20 


22 |. At the same time, several experiments suggested that in some cases the 


establishment of morphogen gradients is associated with nonlinear degradation mechanisms 
when thepresence of signaling molecules self-enhances or self-cat alyzes its removal from the 
system 30|-l33|. Theoretical investigations of temporal evolution of the morphogen gradi¬ 
ents with nonlinear degradation suggested that in this case the local accumulation times, 
in contrast to linear degradation, scale quadraticaly with the distance from the source [30 1. 
But the presented mathematical anal ysis was rather very complicated, and only bounds for 
LAT in several cases where obtained 30 |. 


These observations raised several interesting and important questions concerning the role 
of the degradation in regulating the concentration prohles of signaling molecules. Why 
the degradation accelerates the relaxation to the stationary state for linear degradation? 
Why the actions of linear and nonlinear degradation processes are so different? What is the 
physical mechanism of degradation? Recently, one of us proposed an idea that might resolve 


some of these issues 


22| . It was suggested that the degradation acts as an effective potential 


that pushes signaling molecules away from the source region. It means that the degradation 
will make the diffusion of morphogen molecules effectively biased. However, only qualitative 
arguments have been presented. 


In this paper, we extend and generalize the original idea that the removal of signaling 
molecules works as the effective potential. A new quantitative approach that provides a 
microscopic view on the role of degradation in the formation of morphogen gradients is 
developed. It allows us to explain the differences between various degradation mechanisms. 
We argue that the linear degradation corresponds to a strong potential, leading to strongly 
biased motion of the signaling molecules. At the same time, the non-linear degradation 
creates a potential that is too weak to modify the underlying random-walk scaling behavior 
of the system, affecting only the magnitude of fluctuations. 












II. THEORETICAL METHOD 


Let us start the analysis of degradation mechanisms by introducing a discrete SDD model 
as presented in Fig. la. The cells in the embryo are represented as discrete sites n > 0 
on this semi-inhnite lattice. The signaling molecules are produced at the origin (n = 0) 
with a rate Q. Then morphogens diffuse along the lattice with a diffusion constant D. At 
each lattice site n the molecule can be degraded with a rate kn- It is convenient to adopt a 
single-molecule view of the process where the local concentration of signaling molecules is 
proportional to a probability to hnd the morphogen molecule at a given location 22|, |29 |. 


One can dehne then Pn{t) as a probability of hnding the morphogen at the site n at time t. 
These probabilities evolve with time as described by a set of master equations, 
dPnit) 


dt 


— DPn+l{t) + DPn-lit) — {2D -|- kn)Pn{t), 


for n > 0; while at the origin (n = 0) we have 

dPo{t) 


dt 


= Q + DP,{t)-{D + ko)Po{t). 


( 1 ) 


( 2 ) 


The situation when the degradation rate kn is independent of the concentration of sig¬ 
naling molecules corresponds to linear degradation since the total flux that removes mor¬ 
phogens from the system [knPn{t)] is proportional to the concentration. For the case of 
constant kn = k, this discrete SDD model with linear degradation was fully analyzed before 


22l |. In a more general scenario, the degradation rate might depend on the local concen¬ 
tration, kn = kP{l^~^, where a parameter m specihes the degree of non-linearity, and this 
corresponds to non-linear degradation processes. However, it is not feasible generally to 
obtain full analytic solutions for these non-linear degradation models (with m > 1). 

The main idea of our approach is that degradation acts as an effective potential. This 
suggests that the original reaction-diffusion process with degradation is equivalent to a 
biased diffusion process in such potential but without degradation, as shown in Fig. lb. To 
explain the origin of this potential, let us consider the system in the steady-state limit when 
a stationary non-uniform prohle Pn is achieved. The degradation leads to a concentration 
gradient between any two consecutive sites, and this gradient can be associated with a 
difference in the chemical potentials of the morphogens. 


Hn+i - /Xn = ksTliiPj^^ - ksT In pjf\ 


( 3 ) 








This can also be viewed as an effective potential that influences particles that are not de¬ 
graded. It follows then that this potential can be evaluated as 


U:" = k^TlnP:. 


( 4 ) 


The above arguments indicate that dynamics of the reaction-diffusion model (Fig. la) can 
be mapped into the biased-diffusion model (see Fig. lb), which is much simpler to analyze. 
For the equivalent biased diffusion model we dehne n„(t) as the probability of hnding a 
particle at position n at time t. These probabilities are also governed by corresponding 
master equations, 

= rn+illn+i{t) + - (r„ 5f„)n„(t), (5) 


for n > 0; and 


dt 


= Q + rin,(i) - gono(i)> 


( 6 ) 


for n = 0. The diffusion rates Qn and are related to eac 


1 other via the effective potential 


as can be shown using the detailed balance arguments |34| : 


9n 

Tn+l 



( 7 ) 


This is an important result because it directly couples the original SDD model with degra¬ 
dation to the new biased-diffusion model without degradation. 

One more step is needed in order to have comparable dynamic behaviors in both models. 
The average residence times for the particles at each site provide a measure of relevant time 
scales in the system. It seems reasonable to require that these quantities to be the same in 
both models, leading to 

9n~\~ 'f'n = 2/1 -|- kn- ( 8 ) 


Note that Eqs. 0 and [S] uniquely dehne forward and backward rates in the biased-diffusion 
model. 

To understand the mechanisms of formation of the morphogen gradients the relaxation 
dynamics to a stationary-state behavior needs to be investigated. This can be done by 
analyzing the local accumulation times which are dehned as times to reach the stationary 
state concentration at given position n. The general approach for computing LAT is known 
[^ . but analytical results can only be obtained for the linear degradation model (m = 1). 










We propose to use mean first-passage times (MFPT), which are dehned as times to reach 
a given site for the hrst time, as a measure of dynamics of the formation of morphogen 
gradients. It was shown before that MFPT approximate very well LAT at large distances 


from the source, i.e., for large n 


22 


35| . In addition, the hrst-passage analysis provides 


more clear physical view of the underlying phenomena in the development of morphogen 
gradients. 

Thus, our method of evaluating the formation of signaling molecules prohles consists 
of three steps. First, from the original SDD model with degradation the stationary-state 
prohles are obtained, from which the effective potentials are explicitly evaluated. In the 
second step, the transition rates in the equivalent biased-diffusion model without degradation 
are computed. Finally, these rates are utilized for calculating the hrst-passage dynamics in 
the system. It is important to note here that this procedure is not exact since it involves 
several approximations. 


III. LINEAR DEGRADATION 


To test our theoretical approach, we start with the simplest linear degradation model 
where all dynamic properties are analytically calculated for all sets of parameters 20, l22 |. 


The stationary-state prohle for the SDD model can be easily evaluated 


22 |, 


p{s) 

^ n 


2Qx^ 


(9) 


with X = {2D + k — \/k‘^ + AkD)/2D. This expressions allows us to estimate the ehective 
potential due to degradation for the equivalent biased-dihusion model, 

jjDf 

—— ~ n\nx. (10) 

ksT 


This potential is linear with a slope that depends on dihusion and degradation rates. It is 
also shown in Fig. 5. Employing these results in Eqs. 0 and [HI we obtain the following 
expressions for the forward and backward transition rates. 


9n = g = 


2D + k 

X -f 1 ’ 


Tn+l = r = X 


2D + k 

X -f-1 


( 11 ) 


Note that these rates are independent of the position and the production rate Q. 















In the final step, first-passage dynamics can be evaluated by using known expressions for 
MFPT y, 


= 


TiTi-i-.-rj+i 


n—1 i 

EE 

{x + 1) [x^x^ “ 1) ~ ~ 1)] 


( 12 ) 


(2D + k) (x- ly 

It can be easily checked that in the special case of no degradation in the original system, 
k = 0, this formula reduces to ~ v?/2D at large distances, as expected for a simple 
unbiased random walk. 

It is possible to compare the obtained mean first-passage times from Eq. flT^ with 
analytical expressions for LAT and for MFPT in the original SDD model which are available 
22| . But it is more convenient first to do it in two different dynamic regimes. In the case 
when the degradation rate is much faster than diffusion, k ^ D, it can be shown that 
X ~ D/k, which leads to ~ n/k. It is in excellent agreement with exact results for LAT 


and MFPT for the original SDD model in this limit |22|], tn = ~ {n + l)/k. In the 

opposite limit of very fast diffusion {D ^ k), we have a: ~ 1 — \/kjD and Eq. flT^ yields 


~ n/ VkD. 


(13) 


Exact expressions for LAT and MFPT for the original SDD model give us 22| . 


f ~ 


2k 


1 + 


n + 1 


SDD 




~ nj2VDk. 


(14) 


Thus, for large n our method still correctly reproduces the linear scaling in the local accu¬ 
mulation times, but the amplitude deviates in two times. 

The comparison between predicted MFPT for the biased-diffusion model and for LAT of 
the original SDD model for general sets of parameters is given in Fig. 2. One can see that 
our method approximates the dynamics of the formation of morphogen gradient reasonably 
well. The agreement is better for larger degradation rates where the effective potentials are 
stronger. At the same time, for weaker degradation rates there are deviations, although the 
qualitative behavior is correctly captured. This is a remarkable result given how simple is 
the theory and that it involves several strong approximations. This also suggests that the 
method can be reliably applied to more complex systems with non-linear degradation. 













IV. NON-LINEAR DEGRADATION 


Here we apply our method for systems where the formation of signaling molecules prohles 
is accompanied by the non-linear degradation with = kP^~^ for m = 2,3,.... To eval¬ 
uate the effective potential we need to estimate the stationary-state concentration prohles. 
However, it is not possible to calculate them analytically for general non-linear discrete SDD 
models. But we can use the fact that in the continuum limit {D ^ /c„) the original mas¬ 
ter equations ([H) and ([2]) can be written as the corresponding non-linear reaction-diffusion 
equations, 

( 15 ) 




dt dn^ 

with the boundary condition at the origin 

dP 

n=0 Q- 

on 

These equations can be solved in the steady-state limit, producing 

1 


( 16 ) 


p(«) ^ 


(1 + n/X) 


2 ? 
m — 1 


( 17 ) 


where the parameter A is given by 


A = 


m 


(2T))™(m + 1) 


m+1 


(18) 


kQ^-^ 

It can be shown that the continuum Pn'^ describes also quite well the stationary-state be¬ 
havior of the general non-linear discrete SDD models at large distances from the source. 
This allows us to approximate the effective potentials for non-linear degradation as 

2 


knT 


m — 1 


hi(l + n/X). 


( 19 ) 


This potential is logarithmic, and the degree of non-linearity determines its magnitude as 
illustrated in Fig. 5. It is important to note here that these potentials are always weaker 
than the potential for the linear degradation: see Fig. 5. 

Now using Eqs. [7] and [8] one can obtain the expressions for transition rates in the biased- 
diffusion model. 
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( 20 ) 

















In the final step, again utilizing the analytical framework for the first-passage processes 34j . 
we derive the explicit expressions for the mean first-passage times that approximate the 
formation of morphogen gradients with nonlinear degradation, 


Tn = 


n—1 i 




1 

15 


n—1 


j=0 


2 

m— 1 


i=o j=o 9i9i-i---9j+i9j - 

It can be shown that this expression asymptotically at large distance approaches to 

1 ) 


( 21 ) 


(m 


( 22 ) 


" (m -|- 1) 2D 

This is an important result since it predicts a quadratic scaling for all non-linear degradation 
mechanisms with m > 1. Furthermore, as expected, for very large m, which corresponds to 
effectively no degradation, this formula reduces to a simple random walk dependence. 

Our theoretical estimates for the relaxation dynamics in the establishment of the mor¬ 
phogen gradients for various models with non-linear degradation are presented in Fig. 3. 
One can clearly see that the predicted local accumulation times approach the quadratic 
scaling for large n for all possible ranges of diffusion and degradation rates. The approach 
is faster for larger m. The scaling is independent of the degradation mechanisms, and only 
the amplitude is determined by the degree of the non-linearity m. 

We also compared theoretical predictions with numerically exact values of LAT for dif¬ 
ferent non-linear degradation models. The results are presented in Fig. 4. A remarkable 
agreement between predicted and exact relaxation times is found for m = 3. It can be 
seen that increasing the strength of the degradation (larger k) improves the agreement even 
for small distances from the sources. For m = 10 our theory also works qualitatively well, 
although there are bigger quantitative deviations. It correctly describes the scaling, and 
increasing the degradation rate k decreases the magnitude of these deviations. 

Analyzing results given in Figs 3 and 4, we can make several conclusions about the 
applicability of the developed theoretical method for analyzing nonlinear degradation. Our 
approach correctly finds the quadratic scaling in the local accumulation times. It works 
better for large distances because it calculates only the arrival times which are always smaller 
than the correct LAT that also must include some local rearrangements. At large distances 
the contribution from MFPT to LAT becomes dominant 2^. One can also observe that 


our method works better for stronger degradation, which corresponds to small m values 
and/or large degradation rates k. Most probably, this is due to the fact that our approach 






neglects particle fluctuations that are present even in the absence of degradation. For strong 
degradations these fluctuations become less relevant. 


V. SUMMARY AND CONCLUDING REMARKS 

We developed a new theoretical approach to analyze mechanisms of degradation in the 
formation of signaling molecules prohles during the biological development. The method is 
quite simple, and it provides a full analytical description for all ranges of parameters. It is 
based on the idea that degradation is similar to the effective potential imposed to morphogen 
molecules. The potential pushes signaling molecules away from the source region. It allows 
us to map the original reaction-diffusion process into the biased-diffusion model without 
degradation, which is much easier to analyze. Finally, utilizing the hrst-passage approach, 
the dynamics of relaxation to stationary morphogen gradients can be fully described. 

Despite the fact that our approach involves several strong approximations, it works re¬ 
markably well for different models with degradation. We correctly predict the scaling behav¬ 
ior for the local accumulation times in all cases. As we found for both linear and non-linear 
degradation processes, theoretical method is almost exact for large distances from the source 
and for faster degradation rates. At the same time, for close distances and for slower degra¬ 
dation rates the agreement is mostly qualitative, although the deviations are relatively small. 
The effect of the distance can be explained by recalling that in our method hrst arrival times 
are computed. The correct LAT involve local rearrangements which become less important 
for large distances. The strength and the speed of degradation influence our results because 
the theoretical method neglects the local particle fluctuations due to underlying random walk 
dynamics. These fluctuations are expected to contribute signihcantly to dynamic properties 
for weak and slow degradations, while they are much less important for strong and fast 
degradations. 

The advantage of our method is not only the fact that it gives a fully analytical descrip¬ 
tion of the complex processes during the development of the morphogen gradients. It also 
provides clear physical explanations for the observed phenomena. We can understand now 
why linear and non-linear degradation lead to very different dynamic behaviors. For linear 
degradation we predict that the effective potential is very strong (Fig. 5). The morphogens 
are strongly pushed away from the source region, and as a result a driven diffusion with the 


expected linear scaling is observed. For non-linear degradation processes the effective po¬ 
tentials are much weaker (logarithmic versus linear — see Fig. 5). The particles are moved 
preferentially in the direction away from the source region, but the underlying random-walk 
dynamics is not perturbed much. As a result, the quadratic scaling is predicted and the 
effect of the potential only shows up in the magnitude of fluctuations. These finding also 
suggest that the degradation might be an effective tool for tuning the complex biochemical 
and biophysical processes in biological development. 

Although the presented method captures main features of the degradation processes dur¬ 
ing the formation of morphogen gradients, it is important to note that our approach is 
oversimplified and it involves many approximations. It will be important to test the pro¬ 
posed ideas with more advanced theoretical methods as well as in the extensive experimental 
studies. 
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Fig.l. a) Schematic view of a discrete synthesis-diffusion-degradation model with unbiased 
diffusion, b) Schematic view of equivalent biased-diffusion model without degradation and 
with modified diffusion rates. Lattice sites correspond to embryo cells. 

Fig.2. a) Ratio of the calculated mean hrst passage times and the exact analytical results 
from the SDD model with linear degradation as a function of the distance from the source. 
Different curves correspond to different values of the degradation and diffusion rates, b) The 
same ratio as a function of the ratio of the degradation rate over diffusion. Distance from 
the source is set to n = lO'^, which exceeds the decay lengths for all values of degradation 
rates. 

Fig.3. Theoretically calculated mean hrst passage times as a function of the distance from 
the source for different degrees of non-linearity and for different values of the degradation 
rates: a) m = 2; b) m = 10. 

Fig.4. Theoretically calculated mean hrst passage times (solid circles) and the numerical 
exact results from the SDD model (open circles) for the local accumulation times as a 
function of the distance from the source for diherent degrees of non-linearity and for diherent 
values of degradation rates. For all calculations D = 1 is assumed. 

Fig.5. Ehective potentials acting on morphogens due to degradation. Linear degradation 
corresponds to m = 1, while m = 3 and m = 10 describe diherent cases of non-linear 
degradation. For all calculations k = D = Q = 1 was assumed. 



(a) Synthesis-diffusion-degradation model 



(b) Biased-diffusion model 
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(a) m=2 (b) m=10 
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